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This paper discusses and elaborates on work done in collaboration with 
Drs. T. Stanev of the Bartol Research Foundation and A. K. Harding and J. 
J. Barnard of Goddard Space Flight Center regarding the physics and 
astrophysics of high energy neutrino production and interactions in close 
X-ray binary systems. These studies were stimulated by recent 
observations of ultrahigh energy Y-rays and possibly other ultrahigh 
energy particles coming from the directions of Cygnus X-3 and other 
binary systems and possessing the periodicity characteristics of these 
systems. These experimental reports have potentially profound 
Implications for particle physics, cosmic-ray physics and astrophysics. 
The implications are that such systems are copious sources of high-energy 
cosmic rays, the first such sources discovered. In addition, only a few 
of these systems could account for a significant fraction of the 
ultrahigh energy cosmic rays in the galaxy. Thus, theoretical and 
observational studies of such systems will enable us to understand at 
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least some of the mechanics involved in generating cosmic rays. In 
addition, reports of underground muons produced by neutral particles 
coming from the direction of Cygnus X-3 have raised the possibility of 
new particle physics. 

The Y-ray observations alone imply the production of high energy cosmic 
rays in such copious quantities as to have possibly profound effects on 
the evolution and stability of the binary systems themselves, at least In 
some cases. This raises an interesting new astrophysics problem. On the 
other hand, a detailed understanding of the physics and astrophysical 
effects of high-energy particles, or beams of such particles, enables us 
to determine upper limits for the generation of cosmic-ray particles In 
such binary systems 1 This, in turn, has important Implications for 
Interpreting the underground muon results, which we will discuss. 

We have considered such systems in which a compact object, such as a 
neutron star, is a strong source of high energy particles which, in turn, 
produce photons, neutrinos and other secondary particles by interactions 
in the atmosphere of the companion star. The highest energy neutrinos 
are absorbed deep in the companion and the associated energy deposition 
may be large enough to affect its structure or lead to its ultimate 
disruption. We have evaluated this neutrino heating, starting with a 
detailed numerical calculation of the hadronic cascade induced in the 
atmosphere of the companion star. For some theoretical models, the 
resulting energy deposition from neutrino absorption may be so great as 
to disrupt the companion star over an astronomically small timescale of 
the order of 10,000 years. Even if the energy deposition is smaller, it 
may still be high enough to alter the system substantially, perhaps 
leading to quenching of high-energy signals from the source. Given the 
cosmic-ray luminosities required to produce the observed y-rays from 
Cygnus X-3 and LMC X-4, such a situation may occur in these sources. 

Variability may occur on a much smaller timescale, owing to a strong 
stellar wind generated in the atmosphere of the companion star by cosmic- 
ray energy deposition. Such variability may also be produced by the 
effects of Roche lobe overflow generated by the stellar wind and by 
changes in the internal structure and stability of the companion star 
owing to the absorption of high energy neutrinos. 



3 


It should be clear from all of these considerations, that theoretical 
studies such as the ones discussed here are important for planning muon 
and neutrino astronomy observations with detectors such as DUMAND and 
those detectors at Gran Sasso and Lake Baikal as well as for interpreting 
the results of these observations. 

Reports of air showers with E > 10^ eV from Cygnus X-3 2 , LMC X-4 3 , Vela 
X-l*% and Hercules X-l 5 , have been interpreted as requiring production of 
neutral secondaries by cosmic rays accelerated by the compact partner of 
these systems. If neutral pions are the source of photons that produce 
the observed air showers, then charged pions must also be produced, and 
they will give rise to neutrinos 6 . There is a question as to whether the 
muon content of these showers is high, indicative of hadronically induced 
showers 7 , or low, indicative of electromagnetically (photon) induced 
showers (Klfune, et a!. . Ref. 2). It has been proposed 6 that the 
atmosphere of the companion star provides the grammage required to stop 
the accelerated particles and produce high energy y-rays. (The X-ray 
light curve Indicates a much smaller grammage in circumstellar gas.) 
Reports of underground muons from secondary ultrahigh energy particles 
produced by Cygnus X-3 9 although in question 16 have led to speculations 
of exotic particles and production mechanisms 11 (see, however Ref. 12) 
which would also be accompanied by charged pion-neutrino production, so 
that neutrino beams would also arise In these scenarios. 

What happens when the companion star acts as a beam dump for its compact 
partner which is a powerful primary cosmic-ray accelerator? What is the 
resulting neutrino flux and how much energy is absorbed? How does the 
system respond to rates of energy deposition which may be much greater 
than the stellar luminosity of the companion? These questions can arise 
even independently of the y-ray production model, since it is hard to 
imagine a situation within a close binary system where cosmic rays 
generated at the compact object will avoid hitting the companion star. 
They have been the subjects of our studies. 1 We reconsider them here. 

It is obvious that the "particle beam" produces a hadronic cascade In the 
outer layers of the companion This cascade ultimately gives rise to 
neutrinos and photons. Using previously developed techniques 13 , we can 
calculate the rate of energy deposition in the companion by the cascade. 
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with emphasis on the neutrinos which penetrate deeply into the star. We 
finally conclude with a discussion of the consequences of this energy 
deposition for the stability of the system. 

Hlllas^ has calculated the electromagnetic cascade produced by collisions 
of the accelerated protons in a distribution of gas around the 
companion. He estimated that the observed air shower signal from Cygnus 
X-3 requires a cosmic-ray luminosity for the source of 


L cr - 3 x 10 39 (Ml)M (d /i2 kpc) 2 ergs s" 1 (1) 

Y 

where D y is the duty cycle for production of the photon signal, an Is the 
solid angle into which the high energy particles are beamed and d Is the 
distance from the observer to the source. The power for other sources 
may be similarly estimated from their observed signals, and these results 


TABLE 1 


Aspect Ratio (a/R) and Estimated Cosmic Ray Luminosity 
for Several Binary Systems 


System 

(d) 

binary 

period 

(days) 

(d) 

companion 
mass (M e ) 

(a) 

a/R 

(d) 

distance 

(kpc) 

(b) 

L cr 

(erg/s) 

Cygnus X-3 

0.19 

<4(c) 

-1.4 

>12 

10 39 

Vela X-l 

8.965 

23 

12 

1.4 

-10 37 

LMC X-4 

1.408 

-19 

3.5 

50 

-10 41 

Her X-l 

1.7 

2.4 

6 

4 

-10 38 


(a) Ratio of separation to companion star radius, estimated from Kepler's 
third law assuming the compact object to have mass 1.4 M 0 . Here R is the 
radius of a zero-age main sequence star of the given mass. 

(b) L cr is estimated as in Eq. (1). 

(c) Assuming the compact object to be a neutron star from recent evidence 
that it is a pulsar 14 , and the companion star radius to be less than the 
orbital radius. For stability, the companion star must not exceed its 
Roche lobe by much, giving a limit of -0.5 M for a main sequence star. 

(d) From Ref. 15 and references therein, except for Cygnus X-3. 
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are summarized In Table 1. An important point is that a photon spectrum 
o 

with dn /dE « E , as observed, can be produced by a parent spectrum that 

T 9 O 

is harder than E , including a 6-function spectrum. The E photon 
spectrum is then achieved by electromagnetic cascading, either In strong 
ambient magnetic fields or in the gas. 

In order to understand qualitative aspects of the results it Is helpful 

to see where in the companion star various processes occur. The hadronic 

cascade starts when X < m^/o - 30-40 g/cm in hydrogen for nucleon 

energies around 10^ GeV 1 ^. By ten interaction lengths it has 

saturated. This range of depths occurs in the outer part of the star at 

-7 ^ 

a density of about 2x10 g/cm . In contrast, the interaction length of 
multi-TeV neutrinos is of order 10** g/cm 2 , which occurs deep in the star 
at densities of order 1-10 g/cm . Thus, neutrino production and 
absorption are quite separate. Qualitative estimates as well as detailed 
calculations suggest the following scenario: The hadronic cascade in the 

outer layers of the companion heats and expands the surface and drives a 
stellar wind. These effects are discussed and estimated by Gaisser, et 
al . (Ref. 1). Neutrinos penetrate into or through the star, depending on 
their direction and energy. For close binaries (a/R small) energy 
deposition can be significant. 

We have carried out detailed calculations of neutrino energy deposition 
for two specific companion masses as a function of a/R (Gaisser, et al . , 
Ref. 1), with the results summarized in Table 2. Each case has been 
computed for two extreme spectra, viz. , b = 0 and b = R where b is the 
impact parameter of the trajectory from the compact object through the 
companion (See Fig. 1). In the b * R case, cascading takes place at low 
densities along the edge of the star so that all of the plons and kaons 
which are produced will decay. This gives the hardest neutrino 
spectrum. The b = 0 case, on the other hand, gives the softest neutrino 
spectrum, because the pions cascade to produce lower energy pions which 


decay. These two extremes should bracket the true result, which is 
uncertain because of the extremely violent conditions near the stellar 
surface. 



log E v (TeV) 


Figure 1. Neutrino spectrum at production, given a monoenergetic primary 
cosmic-ray spectrum of 10* 7 eV energy. 
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The dominant process for interaction of the v^'s produced when charged 
pions decay is the charged current mode 


v + N -> u + hadrons. 
p 


( 2 ) 


The cross section for this process rises linearly with neutrino energy 
until the 10-100 TeV range, after which it Increases logarithmically 17 . 

ii o 

The mean thickness of a 2.8 M Q companion star is - 2 x 10 g/cm 
which is equal to the neutrino interaction length for neutrino energy - 2 
TeV. Most neutrinos above this energy will be absorbed. Since the 
critical energy below which neutrinos will be produced by plon decay in 
the outer part of the star Is » 2 TeV, we expect absorption of 
neutrino energy deep in the star to be important. Because of the high 
density of the region in which most neutrinos interact, their energy 
deposition will be fairly localized relative to the point of absorption, 
even for that part of the energy that goes into muons. Charged pions 
will interact and dissipate energy by cascading rather than decaying. 

High energy neutral pions will decay to high energy photons, which will 
produce electromagnetic cascades. Thus, virtually all of the energy of 
the neutrinos which interact will be in the form of low energy photons 
heating the interior of the star. 

Using the standard stellar atmosphere profiles and taking account of pion 
cascading, we have numerically calculated the neutrino energy 
deposition, n = fL , as a function of a/R for two examples, the 2.8 M 

V V 0 

example referred to above and a main sequence star of 15 M . Table 2 
displays the results as a fraction of the cosmic-ray power that is 
deposited in the companion by high energy neutrinos for the two extreme 
production spectra calculated at b = 0 (soft) and at b * R (hard). Here 
L v is the total produced neutrino luminosity and f is the fraction of the 
neutrino power that is absorbed in the companion. Thus, the quantity 
tabulated is n v /L cr . 
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(15 M ) we show the b =■ R production spectrum. The primary cosmic-ray 

0 fl 

spectrum was assumed to be a 6-function with E = 10 GeV. 


TABLE 2 

Neutrino Energy Deposition in the Companion Star as a Fraction 
of Total Cosmic Ray Luminosity 


a/R 2.8 W 15 M Q 



b«0 

b*R 

b=0 

b*R 

1.1 

.0078 

.040 

.011 

.058 

1.2 

.0064 

.033 

.0088 

.046 

1.6 

.0035 

.017 

.0047 

.024 

2.0 

.0022 

.011 

.0035 

.015 

3.0 

.00098 

.0047 

.0013 

.0063 


Let us consider the dynamical effects of the resulting neutrino 
heating. There are at least two relevant quantities to which the 
neutrino heating rate n v should be compared, viz., the Eddington 
luminosity, L^, and the stellar luminosity, L, due to nucleosynthesis. 

If the neutrino energy deposition rate exceeds Lf the whole system may be 
disrupted. If n v Is less than L f but greater than L the system may be 
modified substantially In a time shorter than the lifetime of the system 
as estimated from the observed orbital period derivative. Stecker, et 
a!. 1 , have suggested that the companion may expand In response to the 
neutrino heating and quench the accelerator. We now discuss these 
possibilities in greater detail, in light of the more accurate estimates 
of neutrino energy deposition given in the preceding section. 

If the companion star can absorb a significant fraction of its binding 
energy, it will alter its structure. The time to absorb a binding 
energy. Eg * GM 2 /R, at a rate equal to the heating by neutrinos is 
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at which point the star will be completely disrupted. For heating rates 
near the Eddington limit, this time is relatively short. 

If the neutrino heating rate in the Interior of the star exceeds the 
intrinsic stellar luminosity, the central temperature will initially 
increase. The star will adjust to restore hydrostatic equilibrium, 
reducing the central temperature on a dynamic timescale tp = 1//Gp ~ 1 
hr, which is short compared to tg. In a dynamic timescale, the central 
temperature will increase by only a small fraction of the initial 
temperature. Even though the nuclear reaction rates in the core are a 
sensitive function of temperature, their resulting increase will not be 
large enough to cause a thermal runaway. Instead, the star will expand 
on the much longer timescale tg, maintaining quasi -hydrostatic 
equilibrium. As the star expands the central temperature drops. At the 
same time, the radiative cooling at the surface increases (as a result of 
the larger surface area) and the entire star actually cools as a result 
of the extra heating. From the virial theorem and energy conservation, 
one can show that the initial behavior of the radius of the star as a 
function of time during which continuous heating n v > L occurs is 


R(t) . R 0 / [l-t/t 0 ] t < t Q (4) 

where the constant t Q = (6M 2 /R q )/ 6L, fiL = - L, and R Q is the initial 

radius of the star. This equation assumes that all of the absorbed 
energy goes into work against gravity in expanding the star, which is a 
good approximation until the photon diffusion time becomes comparable to 
t Q . At this point, the star begins to lose some energy through radiative 
cooling from the surface, and the expansion will slow down (unless 
> l-£ ). As the companion star expands it will become somewhat more 
transparent to neutrinos than the initial "undisturbed" star for which 
the absorption calculations were done. Thus, the characteristic time 
would be somewhat longer than the estimates given here, but not 
significantly so for close binaries. Even for heating rates 
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engulf the compact object. We speculate that such effects could lead to 
a mechanism for periodically turning the accelerator on and off and which 
might lead to various observable effects in addition to modulation of the 
high energy signal, such as pulsations, radio outbursts. X-ray outbursts, 
etc. 


TABLE 3 

Neutrino Heating vs. Stellar Luminosities 




% 

L 

l e 

t 0 (* r > 


Cyg X-3^ a) 

10 37 

10 35 

3xl0 38 

4xl0 4 


Her X-l^ 

~10 35 

10 35 

3xl0 38 

— 


LMC X-4^ b ) 

3xl0 38 

10 38 

10 39 

6xl0 3 


Vela X-l< b > 

<10 34 

10 38 

10 39 

— 

(a) 

2.8 solar mass model 




0>) 

15 solar mass 

model 




For 

Hercules X-l 

and possibly Cygnus X-3, the 

2.8 M_ example is the 
0 


relevant one, whereas for LMC X-4 and Vela X-l the companion mass is 
somewhat above 15 M . In the first case Lr - 3xl0 38 erg/s and L - 10 3 ^ 

® OQ OO 

erg/s and in the second the corresponding numbers are 10 erg/s and 10 JO 
erg/s. We now use information from Tables 1 and 2 to make the comparison 
with the neutrino energy deposition. For Vela X-l the separation is so 
great and the companion so massive that n v « L. For Hercules X-l, n v - 
L. The other two cases are the most interesting because in both cases 
the neutrino energy deposition is substantially greater than the stellar 
luminosity. For the 2.8 M case, we estimate n * 10 37 erg/s for the 

9 V y. 

neutrino heating rate and the expansion timescale t 0 - t B = 3.8 x l(r 
yr. From equation (4), it follows that the time for the star to expand 
to the radius of the orbit when a/R Q = 1.2 is t = 6.3 x 10 3 yr. In the 
15 M (LMC X-4) case, we estimate a heating rate of » * 3 x 10 38 erg/s, 

0 V 
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giving t Q = 6.2 x 10 3 yr and t = 4.4 x 10 3 yr for a/R Q = 3.5. Table 3 
shows a summary of luminosities and neutrino energy deposition rate for 
the systems which we have considered. For the Cygnus X-3 and LMC X-4 
systems in particular, our calculated neutrino heating rate exceeds the 
intrinsic stellar luminosity, changing the Internal structure of the 
companion star. We estimate that these systems will evolve or change 
their structure on the characteristic timescale t Q . 

Several possible scenarios may follow the neutrino-heating induced 
expansion of the companion star. If the acceleration mechanism is 
quenched when the compact object is engulfed in the atmosphere of the 
companion, then neutrino production and heating will cease, and the star 
will begin to contract. Once the companion has contracted within the 
radius of the orbit, the acceleration might resume, and some kind of 
periodic behavior may characteristics on timescales of 10 3 - 10 4 yr. 
result. The period would be no longer than t 0 , the timescale on which 
stellar radius changes occur, which is probably too long to be observable 
in the Cygnus X-3 and LMC X-4 systems. (The lower limit on the period is 
- tjj). But if this periodic turning on and off of the acceleration 
adjusts itself so that the time-averaged heating of the star is 
approximately equal to L, then the fraction of time during which 

O 

acceleration is on will be f Qn = L/o v -10 for the Cygnus X-3 case. 

The interval of time during which the source is turned on in a period t Q 
is then f on t 0 < 100 yr. If the oscillation period is between t D and 
several years the v-ray observations would indicate an average luminosity 

30 -1 

of the star to be a few percent of -10 ergs s . This is larger than 
what is expected for an intermediate mass main sequence star. 

If the accelerator does not turn off inside the envelope of the companion 
star, then the heating and expansion of the star will continue. On a 
timescale t 0 , the stellar material surrounding the compact object will 
become optically thick, first to the X-rays (depth of - 1 g/cm 2 ) and then 
to the y-rays (depth of - 50 g/cnr). This timescale may be shortened by 
the orbital decay due to frictional dissipation in the stellar 
atmosphere 1 8 . 

Should the compact object then fall to the center of the companion star 
and somehow continue to accelerate high-energy particles, it is amusing 
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to speculate that a “hidden neutrino source" (e.g. Berezinsky*) might 
result, but here again the stellar atmosphere would continue to expand 
and be blown off. Needless to say, in this "swallowed star" 
configuration, no significant y-ray flux or time-variable neutrino flux 
would be expected, so that this case does not refer to any of the real 
observed systems which are the main concern of this paper. 



Figure 4. Lifetime of a main-sequence companion star against disruption 

by neutrinp heating in the Cygnus X-3 system, scaled by the incident 

39 -1 

proton luminosity L^g = L cr /10 erg s , as a function of its mass from 
Eq (16). The upper and lower curves were determined for the two extreme 
neutrino spectra produced at impact parameters b = 0 and b * R. Values 
of .20 and .05 were used for the neutrino production efficiency e in the 
b = R and b = 0 cases respectively. The mass-radius relation for main 
sequence stars is that given in Ref. 21. 
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Since It lies beyond a very dusty region of the Galaxy, no optical 
emission has been seen from the companion in the Cygnus X-3 system. 

Thus, its mass Is very uncertain, in contrast to the other UHE binaries 
whose companions have been detected optically. Requiring the companion 
to be a main sequence star with a radius less than the orbital radius 
gives a mass upper limit of 4 M e (Vestrand and Eichler 6 ). Requiring the 
star to fit inside its own Roche lobe gives an upper limit of 0.5 M 0 for 
a main sequence star 1 9 and 4 M 0 for a He star 20 . Although the detailed 
results presented above for Cygnus X-3 were obtained for a 2.8 M Q 
companion, the time scales will not be qualitatively different for a 
companion mass as low as 0.5 M Q . The characteristic expansion timescale 
is plotted as a function of stellar mass in Figure 4 for the two extreme 
neutrino production spectra. 1 The break around 1.3 M 0 results from a 
corresponding break in the particular mass-radius relation which we have 
used 21 . The lifetime of a 0.5 M 0 star is between - 10 4 yr and - 10 5 yr, 
with the lower value applicable to a star with a dense wind and the 
higher value being an upper limit. This lifetime is a very small 
fraction of the main sequence lifetime of a 0.5 M 0 star and is not 
significantly different from that derived for the 2.8 M 0 star. 

Several factors may influence the above constraints imposed by neutrino 
heating of the companion star. We have assumed that the accelerated 
particles are emitted isotropically. If the cosmic rays from the source 
are sufficiently anisotropic, then neutrinos would not heat the entire 
volume of the companion star. While there are no direct constraints on 
the solid angle of the particle beam in any of the UHE y-ray sources, in 
the case of Cygnus X-3, constraints on the mass loss from the companion 
star may indicate that the beam heats a small fraction of the stellar 
surface area (Stecker, et al. . Ref 1). However, as long as the total 
energy incident on the star remains the same, we do not expect, given the 
near uniform energy deposition over the star (cf. Figs. 2 and 3), that 
the neutrino absorption will change significantly. There is also 
evidence that Cygnus X-3 is a highly variable source at high energies on 
a time scale of « 100 yr. 22 , which would ease the energy constraints on 
the system, since our limits refer to time averaged energy constraints 
over time scale of 100 years or more. 
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Finally, we note that In the case of Cygnus X-3, n v could be as much as 
1.5 orders of magnitude larger without completely disrupting the 
system. For LMC X-4 It could be a factor of three or so larger. Thus, 
the neutrino-induced signal could be as large as 10-30/1000 m^/year, 
rather than the more conservative estimate of ~ one per year In 1000 
m^. This would require 0 y in Eq. (1) to be much smaller than 1/50. 
However, this duty factor for production of photons is very model 
dependent because of the narrow range of thicknesses for which high 
energy photons can be produced without being reabsorbed. A flux as large 
as this limiting value would still not be large enough to explain the 
reported underground signals from the direction of Cygnus X-3 observed 
with exposures of - 10 m yr 10 , however it could be detectable in the 
large area detectors proposed for Gran Sasso 23 and certainly in a 

C O 

detector with an area as large as 10 m as proposed for DUMAND 24 

The implications of our calculations for the underground muon results and 
their explanation in terms of new particles (cygnets) are potentially 
quite interesting. Cygnets (e.g. quark nuggets, photinos, etc.) are most 
likely generated by the primary particle beam with an efficiency e cy g « 
e y . Defining e C yg/e Y = 5* we note that in order to explain the 
underground muon results, we require a cosmic-ray beam power which is a 
factor higher than those used in the calculations above. In the case 
of Cygnus X-3, however, the underground muon data, when taken at face 
value, indicate that the TeV photon and y fluxes are comparable. We may 
therefore be faced with a contradiction if the muon data are taken at 
face value. 
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